A new barrier penetration formula and its application to a-decay half-lives 
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Starting from the WKB approximation, a new barrier penetration formula is proposed for poten- 
tial barriers containing a long-range Coulomb interaction. This formula is especially proper for the 
barrier penetration with penetration energy much lower than the Coulomb barrier. The penetra- 
bilities calculated from the new formula agree well with the results from the WKB method. As a 
first attempt, this new formula is used to evaluate a decay half-lives of atomic nuclei and a good 
agreement with the experiment is obtained. 
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I. INTRODUCTION 

As a common quantum phenomenon, the tunneling 
through a potential barrier plays a very important role in 
the microscopic world and has been studied extensively 
since the birth of quantum mechanics. One of the earli- 
est applications of quantum tunneling is the explanation 
of a decays in atomic nuclei. The quantum tunneling 
effect governs also many other nuclear processes such as 
fission and fusion. In particular, a lot of new features are 
revealed in sub-barrier fusion reactions which are closely 
connected with the tunneling phenomena [l|, 0, H, Q . 

For most of the potential barriers, the penetrability 
can not be calculated analytically @. Among those po- 
tentials for which analytical solutions can be obtained, 
the parabolic potential @, 0] is the mostly used in the 
study of nuclear fusion. By approximating the Coulomb 
barrier to a parabola, Wong derived an analytic ex- 
pression for the fusion cross section which is widely 
adopted today in the study of heavy ion reactions (see, 
e.g., recent Refs. @, [T3|). The parabolic approximation 
works remarkably well both for the penetrability and for 
the fusion cross section at energies around or above the 
Coulomb barrier [Til ] . 

Apparently the parabolic approximation breaks down 
at energies much smaller than the barrier height due to 
the long-range Coulomb interaction. One may calculate 
the penetration probability numerically by using the path 
integral method or the WKB approximation. However, 
it is highly desirable to have an analytical expression for 
the barrier penetrability when one introduces an energy- 
dependent one-dimensional potential barrier [l2j or bar- 
rier distribution functions [H, [3, [H, [H, E3] • 

In the present work, we derived a new barrier pene- 



tration formula based on the WKB approximation. The 
influence of the long Coulomb tail in the barrier potential 
is taken into accout properly. Therefore this formula is 
especially applicable to the barrier penetration with pen- 
etration energy much lower than the Coulomb barrier. 

As a first attempt and a test study, we apply this new 
formula to evaluate a decay half-lives of atomic nuclei. 
For the a decay, the penetrability is usually calculated 
with the WKB approach [H, [H, i m other words, inte- 
grating numerically the wave number within two turning 
points at which the interaction potential is equal to the 
Q- value of the a decay. Wc will show that the present an- 
alytical formula reproduces the experimental results very 
well, especially for spherical nuclei. 

The paper is organized as follows. In Scc.[Tl]we present 
the new barrier penetration formula. The validity of the 
new formula is investigated and its application to a de- 
cays are given in Sec. IIIII Finally in Sec. IIVI wc summa- 
rize our work. In the Appendix, the detailed derivation 
of the new penetration formula is given. 



II. FORMALISM 

When the penetration energy is well below the 
Coulomb barrier, the barrier penetrability formula de- 
rived from the WKB approximation reads, 



P(E) = cxp 



-2 



flout 



^(V(R)-E) dR 



(1) 



where the potential usually consists of three parts, the 
nuclear, the Coulomb, and the centrifugal potentials, 



V(R) = Vn(R) + V c (R) + 



L(L + l) 
2/iR 2 



(2) 
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i?i n and i? ut are the inner and outer turning points de- 
termined by the relation V(R) = E. 
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By approximating V(R) to a parabola with the height 
Vb and the width hw, Eq. ([T]) is reduced as 



P(E) = cxp 



(3) 



which has been widely used in the study of heavy ion 
reactions. 

Because of the long-range Coulomb interaction, the 
Coulomb barrier given in Eq. ([2|) has a long tail and is 
asymmetric. Thus for the penetration well below the bar- 
rier, the parabolic approximation is not valid. We may 
divide the potential barrier into two parts at the barrier 
position R B . The first part of V(R) with R ul < R < R B 
could still be approximated by half of a parabola and we 
need to evaluate the integration in Eq. (fTJ) in the range 
Rb < R < Rout only. For S wave, the integral in Eq. ((T|) 
is evaluated as, 



P(E) =cxp [-(xi+xa)], 



(4) 



with 



2/i 



(V(R) - E) dR 



(Vb - E), 



(5) 



under the parabolic approximation and 

-Hot 



•"'■2 



'Rb 

2kR B 



^(V(R)-E) dR 




ka Vq 



ln[l + e (R °- RB ^ a ] 



(6) 



where k = y/2jIE/h and r = V C {R B )/E. The details of 
the derivation of Eq. ([6]) are given in the Appendix. It 
should be mentioned that in the derivation of Eq. ([B]), a 
Woods-Saxon form is used for V^{R). 



III. RESULTS AND DISCUSSIONS 

In this section, we use the new formula to study the 
typical barrier penetration problem, a decays of atomic 
nuclei. The a decay half-life is related to the decay width 

rb y ussum 



ft In 2 



The decay width T is calculated as [20j 
r = hvSP{Q) = H£P(Q), 



(7) 



(8) 



where v is the assaults frequency of a particle on the 
barrier, S the spectroscopic or preformation factor and 
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FIG. 1: (Color online) The barrier potential between the a 
and the daughter nucleus for 212 Po and 144 Nd. The solid 
curve shows the exact potential V(R) and the dashed curve 
stands for the effective potential given in Eq. (I12|l associated 
with the parabolic approximation Eq. ((5J and the new barrier 
penetration formula Eq. |BJ|. Note that the two curves are 
almost identical to each other. 



P(Q) the penetrability with Q the a decay Q-value. For 
spherical nuclei, £ is parametrized as (20| 



£= (6.1814 + 0.2988v4.- 1/6 ) x 10 19 s" 



(9) 



and the penetrability will be calculated with Eqs. (J4|, 
©, and ©. 

For the a-nuclear interaction, we adopt the Coulomb 
and the Woods-Saxon potentials and parameters pro- 
posed in Ref. [20T ]. 



V C (R) 



2Ze 2 

"J 



3- 



R 2 



R > Rm, 
R < Rm, 



and 



Vn(R) 



V(A, Z,Q) 



1 + exp [(R - R m )/a] ' 



(10) 



111) 



with A and Z the mass and charge numbers of the daugh- 
ter nucleus and Q the a decay energy. The parameters 
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in these potentials and given in Eq. ([9]) were obtained 
by fitting a decay half lives and cross section data for 
several fusion reactions. It can be easily verified that the 
position of the Coulomb barrier Rb is larger than R m 
thus the use of the Coulomb force given in Eq. (|A3|) is 
valid. 



A. Validity of the new formula 

Before the new formula is used to study alpha decays, 
we investigate in details its validity. First we examine 



is also shown for comparison. i?v ~ 9 — 12 fm is the 
radial position outside of which the nuclear part of the 
a-nucleus potential could be neglected (see the Appendix 
for more details). In our calculations, the width of the 
parabolic potential is obtained by fitting the barrier po- 
tential from the inner turning point R\ n to the position of 
the barrier Rb- Unlike the full parabolic approximation, 
the effective potential is asymmetric and coincides with 
the exact potential very well, especially the outer side of 
the barrier which critically influences a decays. 

In order to examine more closely the accuracy of the 
new formula, we list the calculated penetration probabil- 
ities for a decays of polonium isotopes in Table [T] The 
values in the exponential of Eq. ([5]), x\ and x%, calcu- 
lated from the WKB approach, the parabolic approxi- 
mation and the new formula are compared. One finds 
good agreement between the results from the new for- 
mula and the WKB approach. For X2, the average rel- 
ative root mean square deviation is 0.28 %. This tells 
that the present formula could be used with satisfactory 
accuracy in the study the barrier penetration well below 
the Coulomb barrier. 



B. a-decay half-lives 

The new barrier penetration formula is used to calcu- 
late a-decay half- lives of 344 nuclei collected in Ref . [23| . 
The experimental values of a decay half lives are also 
taken from Ref. [H[ except for 215 Po. In Ref. [H[, 
Iog(lJ^/s) = -3.74 for 215 Po while in Refs. 0, [I| [II 

the experimental value is log(T^2 P /s) = —2.75. We take 
the latter value in the present work. The experimental 
values of the a decay half-lives range from 1CP 7 ~ 10 24 s. 



how the effective potential connected with the new for- 
mula Eq. ((6]) is close to the exact one. Two extreme 
examples are chosen for this purpose, 212 Po which has a 
very short half-life 3.02 x 10~ 7 s and 144 Nd which has a 
quite long half-life 7.24 x 10 22 s [Hj]. The barrier poten- 
tial V(R) is shown in Fig. Q] for these two systems. The 
effective potential, 



(12) 



I 

The Q values of the a decays are also taken from Ref. [23[ 
where these values were calculated from the Atomic Mass 
Evaluation by Audi et al. [13] or from the mass table by 
Moller et al. 

The angular momentum L carried by the emitted a 
particle in a ground-state to ground-state a transition of 
even-even nucleus is zero. In odd- A or odd-odd nuclei, 
L could be non zero. Because the information on L is 
absent, in the present work we assume L = for all a 
decays as usually done [H, SO, SH, H, Hil . 

In Fig. the calculated results and experimental val- 
ues for a decay half lives are compared. In order to 
show it clearly, these 344 nuclei are divided into four 
groups, namely, 159 even-even, 72 even-odd (even-Z and 
odd-N), 66 odd-even, and 47 odd-odd nuclei. The ra- 
tios between the calculated and the experimental values 
S a = log 10 (T 1 c / ^ 1 /T 1 E / 2 P ) are presented in Fig. [3 Two 
dashed lines are drawn to guide the eye. One finds that 
most of the calculated results are of the same order of 
magnitude as the experimental values. 

A statistics of the agreement between the calculation 
and the experiment is made and given in Table UT1 Among 
all these 344 nuclei, there are only seven for which the 
calculated a decay half lives deviate by more than two 
orders of magnitude from the corresponding experimen- 
tal values and 93.90% of them agree with experimental 
values within one order of magnitude. 

Our results arc particularly good for even-even nuclei, 
the calculated half lives for 97.48% of 159 even-even nu- 
clei deviates from the experiment by less than one order 
of magnitude. The ratio S a is less than one for 95.83 % 
of 72 even Z and odd N nuclei, 90.91 % of 66 odd-even 
nuclei and 82.98 % of 47 odd-odd nuclei. The angular 
momentum carried by the emitted a particle might not 
be zero for odd- A or odd-odd nuclei. This will bring in 



V eS {R) = { V C (R) + 
Vc(R), 



V B - \^ 2 {R-R B ) 2 . 

Vc(R)-E 
V C (R B )-E 



4Vc(Rb)-E' 



Rb 



< R < R B , 

< R<R V , 



i?v < i? < R 
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FIG. 2: (Color online) Comparison of the calculated (blue 
crosses) and experimental (red dots) values for a decay half 
lives of 159 even-even, 72 even-odd, 66 odd-even, and 47 odd- 
odd nuclei. 
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FIG. 3: (Color online) The ratios between the calculated and 
experimental values for a decay half lives of 159 even-even, 
72 even-odd, 66 odd-even, and 47 odd-odd nuclei. 
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TABLE I: Comparison of the results for the barrier penetra- 
tion probability for a decays in Po isotopes (charge and mass 
numbers of the a emitter are listed in the first and the second 
entries). The meaning of £1,2 is given in Eq. @. The super- 
script "WKB" means the penetrability calculated from the 
WKB approach, "Para" from the parabolic approximation in 
Eq. ©, and "New" from the new formulas Eq. ©. 



z u 

p 


A p 


J a (MeV) 


WKB 

** > 1 


if w = 


xf KB 


Now 


84 


190 


7.64 


4.9808 


5.0816 


34.9523 


35.0751 


84 


191 


7.48 


5.0093 


5.1146 


36.0311 


36.1527 


84 


192 


7.32 


5.0384 


5.1482 


37.1506 


37.2712 


84 


193 


7.10 


5.0896 


5.2031 


38.7753 


38.8944 


84 


194 


7.00 


5.0980 


5.2165 


39.5213 


39.6397 


84 


195 


6.75 


5.1605 


5.2823 


41.5276 


41.6445 


84 


196 


6.66 


5.1664 


5.2931 


42.2564 


42.3725 


84 


197 


6.41 


5.2292 


5.3592 


44.4401 


44.5546 


84 


198 


6.31 


5.2396 


5.3743 


45.3311 


45.4449 


84 


199 


6.08 


5.2957 


5.4338 


47.5227 


47.6352 


84 


200 


5.99 


5.3036 


5.4463 


48.3954 


48.5072 


84 


201 


5.81 


5.3429 


5.4894 


50.2458 


50.3564 


84 


202 


5.70 


5.3585 


5.5093 


51.4094 


51.5193 


84 


203 


5.50 


5.4050 


5.5594 


53.6532 


53.7618 


84 


204 


5.49 


5.3875 


5.5468 


53.7347 


53.8430 


84 


205 


5.32 


5.4245 


5.5875 


55.7496 


55.8569 


84 


206 


5.33 


5.4015 


5.5693 


55.5907 


55.6978 


84 


207 


5.22 


5.4191 


5.5908 


56.9347 


57.0410 


84 


208 


5.22 


5.4007 


5.5769 


56.8997 


57.0058 


84 


210 


5.41 


5.3038 


5.4892 


54.4768 


54.5836 


84 


212 


8.95 


4.1395 


4.3264 


26.3317 


26.4569 


84 


213 


8.54 


4.2585 


4.4498 


28.5307 


28.6533 


84 


214 


7.83 


4.4720 


4.6690 


32.8310 


32.9495 


84 


215 


7.53 


4.5552 


4.7559 


34.8360 


34.9528 


84 


216 


6.91 


4.7391 


4.9445 


39.4763 


39.5896 


84 


218 


6.11 


4.9669 


5.1798 


46.5717 


46.6806 



some errors for these nuclei in our calculation because 
the centrifugal potential is ignored in the present study. 

The deformation influences the a decay life time both 
on the preformation mechanism and on the penetration 
process [2(1 HH, [H, IHj]. In the present work, we have 
assumed the barrier potential to be spherical. In 68 of 
these 344 nuclei, the spherical potential assumption is 
met well (with |/3 2 | < 0.01 for the daughter nucleus [HI). 
In Table UTTl t he calculated and experimental values of the 
a decay half lives for these nuclei are given. The statis- 
tical summary is also shown in the last line of Table [III 
It is found that the new formula gives very good results 
for these spherical nuclei. In most cases, the differences 
between the calculated and the experimental values of 
^og w T 1 / 2 are smaller than 0.5. The root mean square 

deviation between log^Tf^'/s] and \og w \T^J^ /s] is 0.34. 



IV. CONCLUSION 

In the study of barrier penetration in nuclear physics, 
the parabolic approximation is usually adopted because 
an analytical solution exists for the penetrability of a 



TABLE II: A statistics of the ratios between the calculated 
and the experimental values S a = log 10 (T^ fTf 1 ^) for the a 
decay of 344 nuclei. 68 daughter nuclei are spherical (|/32| < 
0.01) and the results for them are given in the last two lines. 



Nuclei 


\s a \ < 1 


1 < |5 a | < 2 


2 < \S a \ < 3 


All 


323 


14 


7 




93.90 % 


4.07 % 


2.03 % 


Even-even 


155 


3 


1 




97.48 % 


1.89 % 


0.63 % 


Even-odd 


69 


2 


1 




95.83 % 


2.78 % 


1.39 % 


Odd-even 


60 


4 


2 




90.91 % 


6.06 % 


3.03 % 


Odd-odd 


39 


5 


3 




82.98 % 


10.64 % 


6.38 % 


Spherical 


68 










100.00 % 


0.00 % 


0.00 % 



TABLE III: Comparison between the calculated and experi- 
mental q decay half lives of 68 nuclei of which the daughter 
nuclei are spherical (with \{h\ < 0.01). The charge and mass 
numbers of the a decay nucleus are listed in the first and the 
second columns. 



z P 


A P 


Q a 


log 10 [T 1/2 /s] 


z P 


A P 


Q a 


log 10 [T 1/2 /s] 






(MeV) 


Cal 


Exp 






(MeV) 


Cal 


Exp 


52 


106 


4.30 


-3.83 


-4.22 


84 


207 


5.22 


7.24 


8.00 


60 


144 


1.91 


23.17 


22.86 


84 


208 


5.22 


7.22 


7.96 


61 


145 


2.32 


17.53 


17.28 


84 


210 


5.41 


6.13 


7.08 


62 


146 


2.53 


15.61 


15.51 


84 


212 


8.95 


-6.58 


-6.52 


02 


148 


1.99 


23.63 


23.34 


84 


213 


8.54 


-5.58 


-5.38 


63 


147 


2.99 


11.32 


10.98 


84 


214 


7.83 


-3.62 


-3.80 


64 


148 


3.27 


9.46 


9.36 


84 


215 


7.53 


-2.71 


-2.75 


04 


150 


2.81 


13.91 


13.75 


84 


216 


6.91 


-0.61 


-0.82 


66 


150 


4.35 


3.09 


3.08 


84 


218 


6.11 


2.56 


2.28 


66 


152 


3.73 


7.09 


6.93 


85 


213 


9.25 


-6.95 


-6.92 


68 


152 


4.94 


1.06 


1.04 


86 


200 


7.05 


0.11 


0.04 


68 


154 


4.28 


4.64 


4.68 


86 


202 


6.78 


1.08 


1.04 


70 


154 


5.47 


-0.33 


-0.36 


86 


203 


6.63 


1.64 


1.83 


72 


156 


6.04 


-1.66 


-1.60 


86 


204 


6.55 


1.94 


2.00 


72 


158 


5.41 


0.89 


0.81 


86 


206 


6.39 


2.56 


2.74 


74 


158 


6.60 


-2.76 


-3.05 


86 


207 


6.25 


3.15 


3.41 


82 


210 


3.79 


16.16 


16.57 


86 


208 


6.26 


3.08 


3.38 


84 


190 


7.64 


-2.51 


-2.62 


86 


214 


9.21 


-6.52 


-6.57 


84 


191 


7.48 


-2.03 


-1.82 


86 


215 


8.84 


-5.63 


-5.64 


84 


192 


7.32 


-1.53 


-1.47 


86 


216 


8.20 


-3.93 


-4.35 


84 


193 


7.10 


-0.80 


-0.59 


86 


217 


7.89 


-3.04 


-3.27 


84 


194 


7.00 


-0.47 


-0.41 


86 


218 


7.26 


-1.02 


-1.46 


84 


195 


6.75 


0.42 


0.79 


87 


215 


9.54 


-6.94 


-7.07 


84 


196 


6.66 


0.74 


0.77 


87 


217 


8.47 


-4.32 


-4.80 


84 


197 


6.41 


1.71 


2.08 


88 


216 


9.53 


-6.59 


-6.74 


84 


198 


6.31 


2.11 


2.26 


88 


218 


8.55 


-4.17 


-4.59 


84 


199 


6.08 


3.08 


3.64 


88 


220 


7.60 


-1.35 


-1.74 


84 


200 


5.99 


3.47 


3.79 


89 


217 


9.83 


-6.93 


-7.16 


84 


201 


5.81 


4.29 


4.76 


89 


219 


8.83 


-4.56 


-4.92 


84 


202 


5.70 


4.80 


5.15 


90 


218 


9.85 


-6.65 


-6.96 


84 


203 


5.50 


5.80 


6.30 


90 


220 


8.95 


-4.51 


-5.01 


84 


204 


5.49 


5.83 


6.28 


91 


219 


10.09 


-6.85 


-7.28 


84 


205 


5.32 


6.72 


7.18 


91 


221 


9.25 


-4.93 


-5.23 


84 


206 


5.33 


6.64 


7.15 


92 


222 


9.50 


-5.20 


-6.00 
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parabola barrier potential. The parabola approxima- wave, 
tion works indeed well both for the penetrability and for 
the fusion cross section at energies around or above the 
Coulomb barrier. But it fails at energies much smaller 
than the barrier height due to the long-range Coulomb 
interaction. 

In the present work, we derived a new barrier penetra- 
tion formula, Eq. ([5]), based on the WKB approximation. 
We took into account the influence of the long Coulomb 
tail in the barrier potential properly. Therefore this for- 
mula is especially applicable to the barrier penetration 
with penetration energy much lower than the Coulomb 
barrier. We have shown that the present analytical for- 
mula reproduces the WKB results very well. 

This new penetration formula is used to calculate a 
decay half-lives of 344 nuclei with the a-nucleus poten- 
tial given in Ref. [2(|. Satisfactory agreement between 
the present calculation and the experiment is achieved. 
For spherical and even-even nuclei, the results are partic- 
ularly good. Therefore, the new formula could be used in 
the study of barrier penetration at energies much smaller x-i 
than the barrier height. Furthermore, we expect that the 
new formula will facilitate the study of the barrier pene- 
trability where one has to introduce an energy-dependent 
one-dimensional potential barrier or a barrier distribu- 
tion function. 



Re 



! / \I^(V N (R) + V C (R)-E) dR 



! I V \I^(V N (R) + V C {R) - E) dR 



flu 
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R v 
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■y/V c (R)-EJl 
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MR) 
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2fi 



2 l R \j^ { y^ R )- E ) dR - 
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It has been verified that when Rv is not very close to 
Rout, \V N (R)/(V C (R) - E)\ < 1, therefore, 



Rv 
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Rb 
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Since the Coulomb potential outside the barrier (R > 
Rb) is well described by [c.f. Eq. (fTU] , 

Vc(R) = (A3) 

the first term in the above equation can be evaluated 
easily as, 



APPENDIX A: DERIVATION OF THE NEW 
PENETRATION FORMULA 



In order to evaluate the integration x 2 in Eq. ([J|, we 
divide the potential between the position of the barrier 
Rb and the outer turning point i? ou t into two parts, Rb < 
R < an d Rv < R < -Rout ■ Rv should be large enough 
so that the nuclear potential vanishes for R > Rv- For S 



= 2 



^V c {R)-EdR 




Vt - 1 



(A4) 



with k = y/2fiE/h and r = Vc(Rb)/E. For the evalua- 
tion of the second term in Eq. (|A2[) . we adopt a Woods- 
Saxon form for the nuclear part of the barrier potential, 



V N {R) 



1 + exp [(R - Ro)/a] 



(A5) 
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and replace y/Vc(R) — E in the denominator by 

y/V C {R B )-E, 

(2) _ /- flv V2M 



flv ^ V N (R) 
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. e (H-flo)/o]j. 



VV^I E 



{R -R B + <zln[l + e (««-«o)/a]} 
hdl + e^ 0- ^ 70 ]. (A6) 



In the above derivation, we have used the fact that 
exp [(Ry — Ro)/a] ^ 1 for a decay and penetration well 
below the Coulomb barrier. Finally, we have an analyti- 
cal expression for x%, 



■r-2 



2kR v 



I — — arcsin \ / — | 



ka Vq 



\n[l + e (R °- RB)/a }. 



(A7) 
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